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Abstract. The main objective of this article is to study the dynamics of the 
stratified rotating Boussinesq equations, which are a basic model in geophys- 
ical fluid dynamics. First, for the case where the Prandtl number is greater 
than one, a complete stability and bifurcation analysis near the first critical 
Rayleigh number is carried out. Second, for the case where the Prandtl num- 
ber is smaller than one, the onset of the Hopf bifurcation near the first critical 
Rayleigh number is established, leading to the existence of nontrivial periodic 
solutions. The analysis is based on a newly developed bifurcation and stability 
theory for nonlinear dynamical systems (both finite and infinite dimensional) 
by two of the authors 1161 . 



1. Introduction 

The phenomena of the atmosphere and ocean are extremely rich in its organiza- 
tion and complexity, and a lot of them cannot be produced by experiments. These 
phenomena involve a broad range of temporal and spatial scales. As we know, both 
the atmospheric and oceanic flows are flows under the rotation of the earth. In fact, 
fast rotation and small aspect ratio are two main characteristics of the large scale 
atmospheric and oceanic flows. The small aspect ratio characteristic leads to the 
primitive equations, and the fast rotation leads to the quasi-geostrophic equations. 
These are fundamental equations in the study of atmospheric and oceanic flows; 
see Ghil and Childress [6j, Lions, Temam and Wang [HI [13], and Pedlosky [23] . 
Furthermore, convection occurs in many regimes of the atmospheric and oceanic 
flows. 

A key problem in the study of climate dynamics and in geophysical fluid dy- 
namics is to understand and predict the periodic, quasi-periodic, aperiodic, and 
fully turbulent characteristics of large-scale atmospheric and oceanic flows. Stabil- 
ity/bifurcation theory enables one to determine how different flow regimes appear 
and disappear as control parameters, such as the Reynolds number, vary. It, there- 
fore, provides one with a powerful tool to explore the theoretical capability in the 
predictability problem. Most studies so far have only considered systems of or- 
dinary differential equations (ODEs) that are obtained by projecting the PDEs 
onto a finite-dimensional solution space, either by finite differencing or by truncat- 
ing a Fourier expansion (see Ghil and Childress [6] and further references there). 

Date: February 7, 2008. 

Key words and phrases. Boussinesq equations, geophysical fluid dynamics, attractor bifurca- 
tion, onset of Hopf bifurcation. 

The work was supported in part by the Office of Naval Research, by the National Science 
Foundation, and by the National Science Foundation of China. 



2 



C. HSIA, T. MA, AND S. WANG 



These were pioneered by Lorenz p~4] [15] , Stommel [25] , and Veronis [27l [28] among 
others, who explored the bifurcation structure of low-order models of atmospheric 
and oceanic flows. More recently, pseudo-arclength continuation methods have 
been applied to atmospheric (Legras and Ghil [11]) and oceanic (Speich et al. [24] 
and Dijkstra [5]) models with increasing horizontal resolution. These numerical 
bifurcation studies have produced so far fairly reliable results for two classes of 
geophysical flows: (i) atmospheric flows in a periodic mid-latitude channel, in the 
presence of bottom topography and a forcing jet; and (ii) oceanic flows in a rect- 
angular mid-latitude basin, subject to wind stress on its upper surface; see among 
others Charney and DeVore [2], Pedlosky [22], Legras and Ghil [11] and Jin and 
Ghil [10] for saddle-node and Hopf bifurcations in the the atmospheric channel, and 
[20] [H [8] [9] Q1J1 [24] for saddle-node, pitchfork or Hopf in the oceanic basin. 

The main objective of this article is to conduct bifurcation and stability analysis 
for the original partial differential equations (PDEs) that govern geophysical flows. 
This approach should allow us to overcome some of the inherent limitations of the 
numerical bifurcation results that dominate the climate dynamics literature up to 
this point, and to capture the essential dynamics of the governing PDE systems. 

The present article addresses the stability and transitions of basic flows for the 
stratified rotating Boussinesq equations. These equations are fundamental equa- 
tions in the geophysical fluid dynamics; see among others Pedlosky [23] . We obtain 
two main results in this article. The first is to conduct a rigorous and complete 
bifurcation and stability analysis near the first eigenvalue of the linearized problem. 
The second is the onset of the Hopf bifurcation, leading to the existence of periodic 
solutions of the model. 

The detailed analysis is carried out in two steps. The first is a detailed study of 
the eigenvalue problem for the linearized problem around the basic state. In com- 
parison to the classical Benard convection problem, the linearized problem here is 
non-selfadjoint, leading to much more complicated spectrum, and more complicated 
dynamics. We derive in particular two critical Rayleigh numbers R Cl and R C2 . Here 
R Cl is the first critical Rayleigh number for the case where the Prandtl number is 
greater than one, and R C2 is the first critical Rayleigh number for the case where 
the Prandtl number is less than one. Moreover, R Cl leads to the onset of the steady 
state bifurcation while R C2 leads to the onset of the Hopf bifurcation. Both param- 
eters are explicitly given in terms of the physical parameters. The crucial issues 
here include 1) a complete understanding of the spectrum, 2) identification of the 
critical Rayleigh numbers, and most importantly 3) the verification of the Principle 
of Exchange of Stabilities near these critical Rayleigh numbers. 

The second step is to conduct a rigorous nonlinear analysis to derive the bifurca- 
tions at both the critical Rayleigh numbers based on the classical Hopf bifurcation 
theory and a newly developed dynamic bifurcation theory by two of the authors. 
This new dynamic bifurcation theory is centered at a new notion of bifurcation, 
called attractor bifurcation for dynamical systems, both finite dimensional and infi- 
nite dimensional, together with new strategies for the Lyapunov-Schmidt reduction 
and the center manifold reduction procedures. The bifurcation theory has been ap- 
plied to various problems from science and engineering, including, in particular, the 
Kuramoto-Sivashinshy equation, the Cahn-Hillard equation, the Ginzburg-Landau 
equation, Reaction-Diffusion equations in Biology and Chemistry, and the Benard 
convection problem, the Taylor problem; see [16[ I17j and the references therein. 
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We remark that the non-self adjointness of the linearized problem gives rises the 
onset of the Hopf Bifurcation. We prove that the Hopf bifurcation appears at the 
Rayleigh number R C2 . As mentioned earlier, the understanding and prediction of 
of the the periodic, quasi-periodic, aperiodic, and fully turbulent characteristics of 
large-scale atmospheric and oceanic flows are key issues in the study of climate 
dynamics and in geophysical fluid dynamics. It is hoped that the study carried out 
in this article will provide some insights in these important issues. 

Also, we would like to mention that rigorous proof of the existence of periodic 
solutions for a fluid system is a normally a very difficult task from the mathematical 
point of view. For instance, with a highly involved analysis, Chen et al. [3J proved 
the existence of a Hopf bifurcation in an idealized Fourier space. 

The paper is organized as follows. Section 2 gives the basic setting of the problem. 
Section 3 states the main results. The proofs of the main results occupies the 
remaining part of the paper: Section 4 recapitulates the essentials of the attractor 
bifurcation theory, Section 5 is on the eigenanalysis, and Section 6 is on the central 
manifold reduction and the completion of the proofs. 

2. Stratified Rotating Boussinesq Equations in Geophysical Fluid 

Dynamics 

The stratified rotating Boussinesq equations are basic equations in the geophys- 
ical fluid dynamics, and their non-dimensional form is given by 

^ = a(AU - Vp) + aRTe - — e x U - (U ■ V)Z7, 
at Ro 



(2.1) 



dT 

= AT + w-(U- V)T, 

ot 



^ divC/ = 0, 

for [x, y, z) in the non-dimensional domain fl = M. 2 x (0,1), where U = (u,v,w) 
is the velocity fields, e = (0, 0, 1) is the unit vector in the z-direction, a is the 
Prandtl number, R is the thermal Rayleigh number, Ro is the Rossby number, T 
is the temperature function and p is the pressure function. We refer the interested 
readers to Pedlosky [53], Lions, Temam and Wang [T3] for the derivation of this 
model and the related parameters. In particular, the term -^e x U represents 
the Coriolis force, the w term in the temperature equation is derived using the 
stratification, and the definition of the Rayleigh number R as follows: 

(2.2) R= 9 -^h\ 

We consider the periodic boundary condition in the x and y directions 

(2.3) (U,T)(x,y,z,t) = (U,T)(x + 2jw/ ai ,y,z,t) 

= (U,T)(x, y + 2fc-7r/a 2 , z, t), 

for any j, k G Z. At the top and bottom boundaries, we impose the free-free 
boundary conditions: 

U1L UV 

(2.4) (T,w)=0, -7^ = 0, t^ = 0- at z = 0, 1. 

oz oz 

It is natural to put the constraint 

(2.5) / udxdydz = I vdxdydz = 0. 

Jn Jn 
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The initial value conditions are given by 

(2.6) (U, T) — (U, f) at t = 0. 
Let 

ff ={(U,T) e L 2 {Q) 4 | div f/ = 0,u; U=o,i= Q,(u,v) satisfies and 1^3]) }. 
ffi ={(J7,T) 6 H 2 {nf C\H\{U,T) satisfies (gSJ) - (T2~5)lk 
ff ={(tf, T) e H \ (u, v, w, T)(-x, -y, z) = (-u, u>, T)(x, z)}, 

if! =i?x n if. 

Let L R = -A-B R : H l ^ ff (resp., #i -> if ) and G : H x -> ff (rcsp., H x -> H) 
be denned by 

Alp = (-PlaAU - — e x [/], - AT), 
Ro 

B R ^={-P[aRTe},-w), 
GW>)=G(1>,il>), 

for any ip = (U,T) € Hi (resp., ffi), where 

G(ih.,ik) = {-P[{Ui ■ V)E/ 2 ], -(C/x • V)T 2 ), 

for any ?/>i = (U\,T\), V2 = {U2 7 T 2 ) G ffi. Here P is the Leray projection to L 2 
fields, and for a detailed account of the function spaces; see among many others 
[26]. 

Remark 2.1. Note that ffi and ff are invariant under the bilinear operator G in 
the sense that 

G(ipi ,ip 2 )eH, for Vi , V>a e ifi ■ 

Hence, ffi and ff are invariant under the operator L R + G. 

Then the Boussinesq equations (I2.1[) - (|2.5p can be written in the following oper- 
ator form 

(2.7) ^=i^ + GW, tP = (U,T). 

3. Main Results 

3.1. Definition of attractor bifurcation. To state the main theorems of this 
article, we proceed with the definition of attractor bifurcation, first introduced by 
T. Ma and S. Wang in [HQ1]. 

Let ff and ffi be two Hilbert spaces, and ffi <^-> ff be a dense and compact 
inclusion. We consider the following nonlinear evolution equations 

(3.1) {^ = L xU + G(u,X), 

[ u(0) = ItQ, 

where u : [0, 00) — > ff is the unknown function, A e M is the system parameter, 
and L\ : ffi — > ff are parameterized linear completely continuous fields depending 
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continuously on A £ K 1 , which satisfy 

— L\ = A + B\ a sectorial operator, 

(3.2) < A : H\ —¥ H a linear homeomorphism, 

B\ : H\ — » H parameterized linear compact operators. 

It is easy to see [7] that L\ generates an analytic semi-group {e tLx } t >o- Then 
we can define fractional power operators (—Lx)^ for any < \x < 1 with domain 
Hp = D((-LxY) such that H^ C H^ if /x a > fx 2 , and H = H. 

Furthermore, we assume that the nonlinear terms G(-,A) : H^ — » H for some 
1 > fj, > are a family of parameterized C r bounded operators (r > 1) continuously 
depending on the parameter A e R 1 , such that 

(3.3) G(u,X) = o(\\u\\ H J, VAel 1 . 

In this paper, we are interested in the sectorial operator — L\ = A + B\ such 
that there exist an eigenvalue sequence {p k } C C 1 and an eigenvector sequence 
{e k ,h k } C Hi of A: 



(3.4) 



' Az k =p k z kl z k =e k + ih k , 

Rep k — > co (A — > oo), 
_ |Impfc/(o + Rep fe )| < c, 



for some a, c > 0, such that {e kj h k } is a basis of -ff. Also we assume that there is 
a constant < 6 < 1 such that 

(3.5) B x : H 8 — > H bounded, VAeK 1 . 



Under conditions (|3.4p and (J375J), the operator —L\ = A+B\ is a sectorial operator. 

Let {S\(t)} t >o be an operator semi-group generated by the equation (|3.1jl . Then 
the solution of (|3.ip can be expressed as ip(t, ipo) = S\(t)ipo, for any t > 0. 

Definition 3.1. ^4 set E C i? is called an invariant set of i3.1\) if S(t)H = E for 
any t > 0. j4n invariant set Y, <Z H of 13.1]) is called an attractor if E is compact, 
and there exists a neighborhood W C H of ' E suc/i i/iai /or any v/jo £ W we have 

Urn distif(-0(i,-0o),S) = 0. 

i — >oo 

Definition 3.2. (1) FFe saj/ £/ia£ i/ie solution of \3.1}) bifurcates from (ip, A) — 
(0, Ao) to an invariant set if there exists a sequence of invariant sets 
{Q\ n } of US. 1\) such that ^ ^a„> linin^oo A„ = Ao, and 

lim max \x\ = 0. 

n— >oo igfij^ 

(2) J/ i/ie invariant sets £l\ are attractors of i3.1\) , then the bifurcation is called 
attractor bifurcation. 

3.2. Main theorems. In this article, we consider two cases: 

(3.6) cr > 1 and R Cl is obtained only at (j, k, I) — (ji, 0, 1), 

(3.7) a < 1 and R C2 is obtained only at (j, k, I) = (j 2 , 0, 1), 
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for some ji , ji £ N, where R Cl and R C2 are defined in (|5.18|) and Ij5.22|) respectively. 
In the above cases, R Cl and R C2 are given by the following formulas: 



R Cl 



j 2 a 2 a 2 Ro 2 j 2 af ' 

2(a + l)(j 2 2 a 2 + ^ 2 ) 3 27r 2 



if a? (a + l)i?o 2 j 2 2 a 2 • 



Remark 3.3. (1) Condition (|3.6[) guarantees that for R « i? Cl , the first eigen- 
value of Lr \h 1 (resp., Lr \g ) is real and of multiplicity two (resp., one); 
see Remark l5.3l 

(2) Condition (|3.7|l guarantees that, for R as R C2 , there exists only one simple 
pair of conjugate complex eigenvalues of Lr \ g crossing the imaginary 
axis; see Lemma 15.61 

(3) Condition (|3.6|) or (|3.7|) can be satisfied easily; see Lemmas 15.41 and 15.51 



Theorem 3.4. Assume H3.6\) . Then the following assertions for Problem H2.1\) - 
i2.5\) defined in H hold true. 

(1) If R < R Cl , the steady state (U,T) = is locally asymptotically stable. 

(2) For R > R Cl , the problem bifurcates from ((U,T),R) = (0,R Cl ) to an 
attractor Hr — S , consisting of only steady state solutions. 




Figure 3.1. Bifurcation from (0, R Cl ) to an attractor for R > R Cl . 



Theorem 3.5. Assume ^3.7\ l and 



Ro 2 < 



a 2 (l + f T)(j 2 a 2 +7r 2 )3- 

The following statements are true. 

(1) For Problem (2.1]) - I[27S\) defined in H, the steady state (U,T) = is locally 
asymptotically stable if R < R C2 . 

(2) For Problem A2.1\) - [2~5)) defined in H, a Hopf bifurcation occurs generically 
when R crosses R C2 . 
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4. Preliminaries 

4.1. Attractor bifurcation theory. Consider (|3.ip satisfying (|3.2p and (|3.3p . 

We start with the Principle of Exchange of Stabilities (PES). Let the eigenvalues 
(counting the multiplicity) of L\ be given by /3i(A), /32(A), • • • . Suppose that 

(< if A < A , 

(4.1) i?e/3 4 (A) < = if A = A , if 1 < i < m, 

I > if A > A , 



(4.2) Ref3j(X ) < 0, if m+l<j. 
Let the eigenspace of L\ at Ao be 

oo 

Ea= U U^' u e ffl I ( La « "/3j(A )) fc w = 0,w = u + iv}. 

l<j<m k—1 

It is known that dimi?o = m. 

Theorem 4.1 (T. Ma and S. Wang [HICLZ]). Assume that the conditions \3.2\i - 
\3.5\) and |7T7P~|7^P hold true, and u = is locally asymptotically stable for h3.1\) 
at A = Ao- Then the following assertions hold true. 

(1) For A > Ao, \3.1\) bifurcates from (u, A) = (0, Ao) to attractors T,\, having 
the same homology as S™ -1 , with m — 1 < dimY*\ < m, which is connected 
if m > 1 ; 

(2) For any u\ 6 u\ can be expressed as 

u\ = v x + o(\\v\\\ Hl ), v\ € E ; 

(3) There is an open set U C H with G U such that the attractor T,\ bifurcated 
from (0, Ao) attracts U\T in H , where Y is the stable manifold ofu — with 
co- dimension m. 

4.2. Center manifold theory. A crucial ingredient for the proof of the main 
theorems using the above attractor bifurcation theorem is an approximation formula 
for center manifold functions; see [15] . 
Let Hi and H be decomposed into 

(4.3) H x = E$ © £ 2 \ H = E$ © E%, 

for A near Ao eK 1 , where E^ are invariant subspaces of L\, such that dimE* < 
oo, E± — E^, E% — closure of E^ in H. In addition, L\ can be decomposed into 
L\ = Ci © £2 such that for any A near Ao, 

I C\ = L\\ E x : E± — > E£, 



(4.4) 



£9 — L\\ p \ : Ek — > £9 



^ -2 - A I E£ ■ ^2 * ^2 > 

where all eigenvalues of L\ possess negative real parts, and the eigenvalues of L\ 
possess nonnegative real parts at A = Ao- Furthermore, with /i < 1 given by (|3 . 3[1 . 
let 

E^ifJ.) = closure of E^ in H^. 

By the classical center manifold theorem (see among others [3 [26]), there exists 
a neighborhood of Ao given by |A — A | < 5 for some 6 > 0, a neighborhood B\ C E$ 
of x = 0, and a C 1 center manifold function <!>(•, A) : B\ — > E^iO), called the center 
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manifold function, depending continuously on A. Then to investigate the dynamic 
bifurcation of ()3.1|) it suffices to consider the finite dimensional system as follows 

(4.5) -£=£ix + gi(x,& x (.x),\), xeB x cE*. 
at 

Hence, an approximation formula for the center manifold function is crucial for 
the bifurcation and stability study. 

Let the nonlinear operator G be in the following form 

(4.6) G(u,X) =G n (u,\) + o(\\u\\ n ), 

for some integer n > 2. Here G n : Hi x • ■ ■ x H\ — > H is a n-multilinear operator, 
and G n (u, A) = G n (u, ■ • ■ , u, A) . 

Theorem 4.2. [16] Under the conditions {j.Sty , \4-4\ ) an d $4-6\) > the center manifold 
function <i>(x, A) can be expressed as 

(4.7) A) = (-C^y 1 P 2 G n (x 1 A) + o(\\x\\ n ) + 0(|Re/3| \\x\\ n ), 

where L\ is as in {4-4\) > ^2 : H — > E 2 the canonical projection, x 6 E± , and 
(3 — (/3i(A), • • • ,/3 m (A)) the eigenvectors of H\. 



5. Eigenvalue Problem 
The eigenvalue problem of the linearized problem of (|2.ip ~ (12.4p is given by 



(5.1) 



a(AU - Vp) + o-RTe ~ ^-e X U = (3U, 
Ro 

AT + w = pr, 
divU = 0, 



supplemented with (|23|) and For ip = (U,T) satisfying flXg]) and IJ23J), we 

expand the field ^ in Fourier series 



(5.2) 



j.k— — oo 



Plugging (|5.2p into (|5.ip . we obtain the following system of ordinary differential 
equations 

o-(Dj k Ujk - ijaipjk) + -jrVjk = /3ujfe, 

a(Dj k Vj k - ika 2 pjk) - = u j*> 

ixo 

(5.3) <( D jk w jk - p' jk + RT jk = o-~ 1 Pw jk , 

DjkTjk + Wjk — @Tjk> 
ijctiUjk + ika 2 v jk + w' jk = 0, 

k Uj k |*=0,1 = Vjk \z=0,l— Wjk |*=0,1 = T.jk |«=0,1= 0) 
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for j, k £ Z, where ' = d/dz, Dj k = d 2 /dz 2 — a 2 k and a 2 k = j 2 a\ + k 2 ot\. If 
Wjk ^ 0, (|5.3p can be reduced to a single equation for Wj k (z): 

(5.4) {(D jk - (3){aD 3k - (5) 2 D ]k 

+ ^( D jk ~ P){D ok + a 2 k ) + aRa 2 jk (aD jk - f3)}w jk = 0, 

(5.5) Wj k = w jk = wf k = wf k ' =0 at z = 0, 1, 

for j, k G Z. Thanks to ()5.5|) . can be expanded in a Fourier sine series 

oo 

(5.6) w jk (z) = ^2 w jki sin lirz, 

i=i 

for (j, k) G Z x Z. Substituting (|5.6[) into (|5.4p . we see that the eigenvalues /3 of 
the problem (|5.1|) satisfy the cubic equations 



/ 2 tt 2 a 2 , 

(5.7) /3 3 + (2a + l) 7 ? H /3 2 + [(a 2 + 2^ + __ - *R-£]f3 

no Ijkl Ijkl 

+ a 2 1 % l -a 2 Ra% + l ^ = 0, 
for j,k G Z and / G N, where 7 2 fci = a 2 fc + Z 2 7r 2 . In the following discussions, we let 

9iu{fi) = 09 + 7|h)[(/3 + ^7i«) 2 + ZVi? - 2 7ifc 2 ], 

(5.8) 09) = aRafaikiiP + ^m), 
fjki{0) = QjkiiP) ~ hjki((3), 

and (3jkii(R), (3jm(R) and /3 jk i 3 (R) be the zeros of /jfcj with 
Reffijm) > Re(J3 Jk a) > Re((3 jm ). 

5.1. Eigenvectors. In the following discussions, we consider the following index 
sets: 

A! = {(j, k, I) G Z 2 x N | j > 0, (j, k) ± (0, 0)}, 
A 2 = {(j, K G Z 2 x {0} | j > 0, (j, fc) ^ (0, 0)}, 

A 3 = {(i,fc,0e{(0 ) o)}xN} 1 
A = Ai U A 2 U A 3 . 

1. For (j, fc,0) G A2, we define 

ipf 3k0 — {ka<2. sin(jaix + ka^y), — ja\ $m(jaix + fea^y), 0, 0)*, 
= (-ka 2 cos(j aia; + ka.2y\joL\ cos(jaia; + ka^y), 0, 0) , 

/?A 2 = u (j,fe,o)eA 2 {/3ifco}, 
where /3jA;0 = ~ a ljko = — <7Q;2 fc = — a (j 2a i + fc 2 a 2 ). It is not hard to see that 
L R {4 ]k °) = Ml™ and L R (^ k °) = ■ 
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2. For (0,0,0 e A 3 , we define 

^oou = (0,0,0,sinZ7rz)*, ip 13 " 012 = (cos Ittz, 0, 0, 0)*, 

^/3oo!3 = (0,coste,0,0)*, E m = spaniip 000 ' 1 ,ifj Po012 ,ifj l3 ° 013 }, 

Pa 3 = Ug x U, =1 {/?ooig}, /?x 3 = Ugi^oa}, 

where /3 oii = ~7ooi = -^ 2 , /?ooi2 = -c7oo; ~ ife* and Am3 = ~^lo i + ife*- :t 
is easy to check that 

L R (^ oon ) =#x>uV /W , 

L R {ip Paal2 ) = -<r-fo m ip f:>om2 - —if, 130013 , 

Ro 

L R (il) f) ° 013 ) = —if; 1300 ' 2 - ajL,^ 0013 . 
Ro 

3. For (j, fc, Z) £ Ai, we define 

= ( — ^~~2 — &ijx(joi.\X + hotly) cos Z7rz, 1 — sin(jaia; + ka^y) cos /7rz, 

cos(jaia; + kotiy) sinZ7rz, 0)*, 



,2 



ka2ln jctxln 
'ikl — i — 9 — sinliaix + ka^y) cos Inz, 5 — sinnaix + fca^y) cos «7T2:, 0, 0), 



jki = (0, 0, 0, cos(jaix + kot2y) sin/7rz) , 



;,4 



; — (- — \ — cos(jaia; + ka 2 y) cosIttz, — | — cos(jaix + ka 2 y) coslnz, 
a jk a jk 



sin(jaia; + ka 2 y) s'mlirz, 0)*, 

g fca^TT jOL\l-K f 

4>iki = ( 9 — cosyaix + ka 2 y) cos f7rz, — 5 — cosuaiir + rea<22/) cos Z7rz, 0, 0) , 

<^ fc ; = (0, 0, 0, sin(jaia; + ka2y) sin Ittz) 1 , 

E )ki = span{(j)] kl ,<j) 2 kl , <j)%i}, E 2 kl = span{4>% u <f> 5 jkl , cj) 6 jkl }, 
Ejki = E] M © /3 Al = U^,^)^ U 9=1 {/3jfei ? }- 

It is easy to check that Ej kl and E 2 kl are invariant subspaces of the linear op- 
erator L R respectively, i.e., L R (Ej kl ) C Ej kl and L R (E 2 kl ) C fyjy. The char- 
acteristic polynomial of (resp., L# \E 2 kl ) 1S given by fj k i as defined in 
(|5.8p . Since _Ej fei (resp.,i?J feZ ) is of dimension three, the (generalized) eigenvectors 
of L R \ EjM , U q=1 {^ k "} (U 3 q=1 {^2 Jk ' q }), form a basis of E) kl (resp., E 2 kl ), i.e., 
span{u3 =1 {^f Jfc!9 }} = 4w (^P., span{lj3 =1 {^ ! * }} = £*,). If jklq is a real 



zero of fjki, the eigenvector corresponding to (3j k i q in (resp., E 2 M ) is given by 

(5-9) ^ k,q = 4>) kl + Ax + 4i(/W$m> 

= <t>% l + A 1 (0 jMq )4> 5 jkl +A 2 ^ jMq )4>% l ), 

where 

(5.10) Ax(/?) = 2 , , A a (/8) = 



i2o(^ + a7? fcJ )' zvw P + 1 2 U 
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If Pjkiq-L = Pjkiq 2 (imaginary numbers) are zeros of fjki, the (generalized) eigenvec- 
tors corresponding to 0jkl qi and (3jkiq 2 in Ej kl (resp., Ej kl ) are given by 



(5.H) 



4 3klqi = <t>) k i + Mf3jki qi )<t>%i + MPjunWjki, 
4 Dklq2 = h (P jk i qi )4> 2 m + h (P jk i qi )$ w , 



4 3klq2 = h{(3 0klqi )^ kl +h{P 3klqi )4>%u 



where 
(5.12) 

The dual vector corresponding to ipi 3k ' 9 (resp., ip2 ,klq ) l& given by 



Ri{(3) = Re(Ai(0)), MP) = Re(A 2 (/3)), 
h(p) = Im(Ai(0)), I 2 (0) = Im(A 2 ((3)). 



(5-13) = 4>) k i + Ci(f3 3k i q )tf kl + C 2 ([3 jklq )<p% u 

{< 3klq = <t>%i + Ci + C 2 (/W#w). 



where 



< 5 -»> = ^ CM -J^,- 

The dual vector ^ Jklq (resp., ^ kl ") satisfies 

(5.15) < Vf jfc ' 9 * , ^>i klq >h= (< Va 3 """* . ^2 kl " >h= 0), 
for g* 7^ g. 

We note that E nklh is orthogonal to £/j 2 fc 2 j 2 for (ji, fci, Zi) ^ (j 2 , k 2 ,l 2 ) and Ej k[ 

is orthogonal to E 2 jkl for (j, fc, l)eAi. Hence the dual vector vf jk,q (resp., ^ kl ") 
satisfies 

(5.16) <^*?' w « >»=0 for ^6(U (j , ifc .,,. mi) £ i . fc . i .)U^ i 
(< > ff = for V G (U( J .iM. )/W )E J . M .)u4 1 ). 

In view of the Fourier expansion, we see that \Jrj tk n^\Ej k i is a basis of H\ and 

(U (3 -,fc,l)€A 1 4« ) U ( U 0'.fe,0)eA 2 {^f 3M }) U (U ( o, ,i)eAs i^ 0011 }) is a basis of ^ . Hence, 
by the discussion above, we have the following conclusions. 

a) The set 0h x = /?Ai U /3a 2 U /3a 3 consists of all eigenvalues of Lr \h 1 , and the 
(generalized) eigenvectors of L R \ Hl form a basis of Hi. 

b) The set 0^ = 0a x U /3a 2 U consists of all eigenvalues of Lb, and the 
(generalized) eigenvectors of Lr \ g form a basis of Hi . 

c) Re{0) < for each & 0a 2 U /3a 3 . 

Lemma 5.1. is small, then Re(0jki q (R)) < /or eac/i 0jhi q S /?Ai- 
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Proof. Plugging = j 2 kl f3* into f jk i, we get f jk i(/3) = l%Jjkl{P*), where 

fMn = (0* + i)G9* + °? + + 1) - vR^iP* + 

Hence, we only need to show that the real part of each zero of fjki is strictly negative 
when R is small. We observe that fjki{(3*) > for all (3* > provided R < 1 + er~ 1 . 
Therefore, if all zeros of fju are real numbers, we are done. 

For the case where only one of the zeros of fjki is real, this real zero, (3*, is a 
perturbation of —1. There exists an e ( depending on a only) such that —(l + 2o~) < 
(31 < provided R < e. This makes the real part of the other two zeros of fjki 
strictly negative and the proof is complete. □ 

5.2. Characterization of Critical Rayleigh Numbers. Based on the above 
discussion, we know that only the eigenvalues in (3^ depend on the Rayleigh number 
R. Hence, to study the Principle of Exchange of Stabilities for problem (|5.1|) . it 
suffices to focus the problem on the set f3\ 1 . We proceed with the following two 
cases. 

Case 1. (3 = is a zero of fjki if and only if the constant term of the polynomial 
fjki is 0. In this case, we have 

_l]u Z 2 7r 2 (a% + n 2 f 7T 2 

(5 ' 17) R a% + a*Ro 2 a% ~ a% + oWafr " 

Hence the critical Rayleigh number R Cl is given by 

<a/ 6 72^2 ~,6 2 

0,fc,0eAi a ife cr ito « ife a jiki a Ho a Jiki 

for some (ji, fci, 1) 6 Ai. 

Case 2. A careful analysis on (|5.7p shows that (3 = ai (a =fi 0), a purely 
imaginary number, is a zero of /jfc; if and only if the following two equations hold 
true: 

1X0 Ijkl Ijkl 

(2^ + l)l%i\^ + 2ahj kl + - *R^} 

1X0 Ijkl Ijkl 



2 6 2 D 2 

= ° Ijkl - a Ra jk 

In this case, we have 

' V fjki 



l 2 TT 2 

Ro 2 ' 



(5.19) R = 



(5.20) R < 



2(^ + 1)7^ 2l 2 n 2 

a% + (a+l)Ro 2 a%> 

(^ + 2)7^ . l 2 n 2 



a 



2 



<jRo 2 c 



Plugging (|5.20p into (|5. 19|) . we derive an upper bound for Ro 2 , 
(5.21) Ro 2 < 
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which could only hold true when a < 1. 

As in Case 1, the minimum of the right hand side of (|5. 19|) is always obtain at 
1 = 1. Hence the critical Rayleigh number R C2 is given by 



■ 2(^+1)7^ , 2Z 2 - 2 



(5.22) R C2 = min { > . ., ,, , ., 

( 3 ,fe,06Ai a 2 jk {a + l)Ro 2 a 2 jk 



2(^ + l)7L 2 i , 2n 2 



< k2 (a + l)Ro 2 al k2 



for some (j*2>&2>1) G Ai. In the case of a < 1, (|5.2ip with I = 1 implies R C2 is 
smaller than R Cl . Hence, for Problem (|2 . 1 [) - f|2 . 5[> . R Cl is the first critical Rayleigh 
number if a > 1 and R C2 is the first critical Rayleigh number if a < 1. Therefore, 
the Principle of Exchange of Stabilities is given by Lemma 15.21 and Lemma 15. 61 



Lemma 5.2. For fixed a > 1 and Ro > 0, suppose that (a 2 k ,l) = (o^fcjjl) 
minimizes the right hand side of \5. 17\ ), then 



< if R< R C1 
(5.23) PnkMfy < =0 if R = R C1 



> if R> R. 



Cl 



(5.24) Ref3 3Hq {R) <Q for (a%,l) ? (a 2 lfel , 1), g = 1, 2, 3, R near R C1 . 

Proof. By the above discussion, we only need to show that the first eigenvalue 
crosses the imaginary axis. We note that fj 1 k 1 i{P) = is equivalent to g J1 fc 1 i(/3) = 
hnkMP), i-e., 

(5.25) (P + ^ kl iW + °7l kl i) 2 + l 2 K 2 fo- 2 7iXi} = ^Ik^XiiP + ^nka)- 



We see that both gj x k x x and hj lkl i are strictly increasing for > — 7 2 lfcll ( since 
a > 1 ). Let Ti be the graph of r\ = gj 1 kii(0) and T 2 be the graph of r\ — hj^iiP) 
as shown in Figure [5Tl When R = R Cl , Point So, the intersecting point of Ti and 
T2 corresponding to /3j 1 fe 1 i(^2) ( i.e., the (3 coordinate of Sq is f3j lkl i(R)), is on the r\ 
axis. When R increases (resp., decreases), So becomes Si ( resp., 5*2). This proves 
(|5.23[) and the proof is complete. □ 
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Figure 5.1. 



Remark 5.3. (1) In the proof of Lemma 15.21 as shown by (|5.25[) and Fig- 
ure [53] we see that, for R w R C1 , the first eigenvalue l3j 1 k 1 n is a simple 
zero of fhki\{0). We have seen in Section 5.1 that there are eigenvectors 

t/>f nfcl " 6 Ejikil and ■02 Jlfcl11 G E? ikil corresponding to fy^u- Therefore, 
the multiplicity of the first eigenvalue of L \h x (resp., L \ ) is mH 1 = 2m 
(resp., = m ), where m is the number of (j,k, l)'s (s Ai) satisfying 
a jk = a iifei- Hence, Condition (|3.6p guarantees that, for i? » i? Cl , the first 
eigenvalue of L# (resp., L^j |g ) is real and of multiplicity two (resp., 
one). 

(2) For the classical Benard problem without rotation, the second term on the 
right hand side of (|5.17j) . hence the second term on the right hand side of 
(|5.18j) . is not presented. Therefore, the first critical Rayleigh number of the 
classical Benard problem depends only on the aspect of ratio; while the first 
critical Rayleigh number of the rotating problem depends on the aspect of 
ratio, the Prandtl number and the Rossby number. And it is clear that 
the first critical Rayleigh number of fast rotating flows is remarkably larger 
than the first critical Rayleigh number of the classical Benard problem. 
This indicates that the rotating flows are much more stable than the non- 
rotating flows. 

(3) R Cl is the first Critical Rayleigh number if the Prandtl number is greater 
than one. For the case where the Prandtl number is smaller than one, R C2 
is the first Critical Rayleigh number and, in general, there are a few critical 
values between R C2 and R Cl ■ 

For x > 0, b > 0, we define 

(5.26) f b{x) = ( - + - 2)3 + & . 

x 

Let x — aj k , then the right hand side of (|5.18[) could be expressed as fb x (x) , where 
bi = a 'i Ro i ; an d the second line of (|5.22[) could be expressed as 2(<r + l)fb 2 (x), 
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where b 2 — /_ , tts d „9 ■ Consider 



((T+1)2_R 2 ' 

(2a;-^ 2 )(a ; + 7r 2 ) 2 



(5.27) / b (x) 

As shown in Figure [5?2l it is easy to see that 

a) for x £ (0, oo), fb(x) has only one critical number Xb, 

b) f b ( x ) < if x < Xb, 

c) f b (x) > if x > x b , 

d) fb(xb) is the global minimum of fb(x), and 

e) Xfc is strictly increasing in b , hence, a;^ > Xb 2 > \- 




Figure 5.2. 



In Lemmas 5.4 and 5.5, we consider the following different conditions 

(5.28) x bl <a\< a%, 

(5.29) a\ < ^x bl < 2x bl < a\, 

(5.30) Xb 2 < ot\ < a 2 , 

(5.31) a\ < —Xb 2 < 2xh, < a\. 

5 

Lemma 5.4. (1) Condition 13. 6\) holds true under the assumption i5.28\) . 
(2) Generically, Condition \3. 6\) holds true under the assumption h5.29j) . 

Proof. (1) Under the assumption (|5.28j) . by c), we conclude that R Cl is only 
obtained at (j,k,l) — (1,0,1), i.e. j\ = 1. 
(2) Under the assumption (|5.29p . there exists j* > 2 such that j* 2 a\ < Xb 1 < 
(j* + \) 2 a\. We note that 

( * + l) 2 a 2 { < 2j * < 2Xbl < " 2 if j * ~ 3 ' 

U ' 1 \ = 9aj < |x 6l < 2x bl < a 2 if f = 2. 
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Hence, by b) and c), we conclude that 

R C1 = mm{f bl (f 2 al)J bl ((f + l) 2 a\)}, 

i.e., ji = j* or ji = j* + 1. Note that, by b) and c), generically / bl 0'* 2 a?) + 
/&i((i* + 1) 2q; i)- The proof is complete. □ 

Lemma 5.5. (1) Condition {3.7jj holds true under the assumption h5.30\ ). 
(2) Generically, Condition \3.ty holds true under the assumption h5.31\) . 

Proof. Consider 

R C2 = min {2(o- + l)/b 2 (a| fc )}- 
The rest part of the proof is the same as the proof of Lemma 15.41 □ 

Lemma 5.6. Assume \3. 7\ ), R « R C2 and Ro 2 satisfies \5.21\) for (j, k, I) — 
( J2 ,0,l), i.e., Ro 2 < J^Z ] S > then{/3 h011 (R),f3 h012 (R)} (f3 h011 (R) = hO12 (R)) 

is the only simple pair of complex eigenvalues of the problem \5.1)) in space Hi sat- 
isfying 

!< if R < R C2 , 
= if R=R C2 , 
> if R > R C2 , 



(5.33) Re[3 jklq {R) < for {a 2 k ,l) + (a% , 1), q = 1, 2, 3, R near R c . 2 . 

Proof. We only need to prove (|5.32|) . Under the assumptions of the lemma together 
with (|5.11|) . (|5.19|) and ()5.20|) . by the discussion in Case (2) at the beginning of this 
subsection, we know that {f3j 201 i(R) , f3j 2012 (R)} is the only simple pair of complex 
eigenvalues of L R \^ with Re(f3j 2 ou(R C2 )) = Re((3 j2 oi 2 (R C2 )) = 0. Since 0j a 013 (R) 
(real), f3j 2 on(R) and (5j 2 o\ 2 (R) are zeros of fj 2 oi, we know that 

P h0 i 3 (R) = -(Re(P h0 u( R )) + Re(f3 h012 (R))) - (2a + 1)7? 01 . 
Hence (15.32)) is equivalent to 



(5.34) 



> -(2u + lMm if R<R, 
h oi 3 ( R ){ = -(2a+ 1)7? 01 if 11 = 11 

hoi 



< -(2(7 + 1)7^ if R> R, 



C2 ! 
C2! 
C2! 



which is true as shown in Figure 13731 This completes the proof. 
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Figure 5.3. 



□ 



Lemma 5.7. For fixed a\, 012 > and a > 1, i? ci — > oo as i?o — > 0. More 
precisely, R Cl = 0(Ro~z). 

Proof. Since b x = ^2^, by (|5"?7j) . z 6l = 0(b\) as i?o -> 0. Hence, 
R Cl = 0(f bl (x bl )) = 0(6?) = O(Eo-t). 

□ 

6. Proof of Main Theorems 

6.1. Center manifold reduction. We are now in a position to reduce equations 
of (|2.1[) - (|2.5[) to the center manifold. For any ip = (U,T) € Hi, we have 

oo 3 

y,fc,i)eAi 9=1 

OO 3 

+ E (»i*o^ M + t&rf* ) + E E aW 00 " • 

(j,fc,0)eA 2 1=1 9=1 

Under the assumption (|3.6|) . the first critical Rayleigh number is given by 

(6.1) R C1 = + 2 . 

a| l0 cr 2 Ro 2 a z jl0 

In this case, the multiplicity of the first eigenvalue is two and the reduced equations 

of (|XT j) -([53 ]l are given by 

(6.2) 
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Here for V>1 = {Ui,T ± ), i> 2 = (U 2 ,T 2 ) and </>3 = (U 3 ,T 3 ), 
G(V>i, ip 2 ) = -{P(Ui • V)(7 2 , (C/i • V)T 2 )< 

and 

"1 r2ir I a2 r2ir/a-L 



[< {Ux'V)U 2 ,U 3 > R3 
+ (Ui ■ X7)T 2 T 3 ]dxdydz, 



where P is the Leray projection to 1? fields. Let the center manifold function be 
denoted by 

(6.3) $= (*i (^011,1/^011)^1 + $2(^1011, Vh0u)^ 2 )' 

/Mft l0 ii 

The direct calculation shows that 

Am 2 



(6.4) 



G(^f Jl011 ,^f Jl011 ) 
G(7/;f Jl011 ,V2 Jl011 ) 



G(^' l011 ,^f Jl011 ) 
G(V>2 Jl01 \^2 Jl011 ) 



2jiai 2 



• cos 2irz. 



-( cos 2nz, 

2jiai 



A l7 r 2 
2jiai 
-A l7 r 2 
2ji«i 



(cos 27T2; — cos 2jiaix), 0, 0)*, 
(cos 2jiaia; + cos27rz), 0, 0)*, 



—Ain 2 . A 2 tt . t 

-(0, sin 2j\a\x, 0, sin 2itz) . 

2ji«i 2 



G(^f Jl011 ,*f 1011 ) 



ftlOll .T.ftjOll 



(6.5) 



G(^' l011 ,*f 1011 ) 



G(V» ; 



> 2 



Gl7T 2 G27T t 

-(0, — sin2jiaix, 0, — — sin27rz) , 

2jiai 2 



) = -( 



2jiai 



COS 27TZ, 



2jiai 



• COS 27TZ 



G l7 r 2 
2jiai 
-G l7 r 2 
2jiai 



) = -(o, 



-G l7 r 2 



(cos27rz — cos2jiaia;),0,0)*, 

(cos 2jia\x + cos 2irz), 0, 0)*, 
G 2 tt 



sin 2jiaix, 0, — — sin 27rz)* , 

2jiai 2 



where = Ai(/3 il0 n), A 2 = A 2 (/3 ilQ n) Gi = Ci{fi jmi ) and G 2 = G 2 (/3,- l0 ii)- 
Hereafter, we make the following convention: 

o(2) = o(4 0U + y 2 l011 ) + 0(| /3 Jl0 ii(-R) | -(4oii + l&oii)). 

o(3) = o((4on + ^ l0 ii) 3/2 ) + 0(| /3 jl0 u(R) I -(4ou + ^on) 3/2 ), 

o(4) = o((4on + ?4oii) 2 ) + 0(| P hQ u{R) | -(x 2 l011 + yf i011 ) 2 ). 

By Theorem |4~2"1 and (f6~4>(|6~l)j) . we obtain 

(6.6) $ = $f< 2 ") TO ^f (23l)0 ° + ^t 2 "' 00 ^ 2 "' 00 + $A>021^0021 + o(2 ) ; 
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where 

$ % 31 ,oo = J^_^ ou -j / | iOU ) + (2), Vf (23l)00 = (0,-2j iai sm2j 1 a 1 z,0,0) t , 

CTa (2jl)0 

$2 (3il) °° = -4^(2^ioii2/ioii) +o(2), V2 (23l)0 ° = (0,2i 1 aicos2aix,0,0) t , 

aa (2j 1 )0 

$ i° 021 = ^(4011+^011) + o(2), V'f 0021 = (0, 0, 0,811127^. 

Note that for any ^, G -ffi(i = 1, 2, 3), 

(6.7) <G(ih,ih),ih> H =0, 

(6.8) < G(ipi,i> 2 ),tl>3 >h= - < G(^i,^3),^2 >h; 
and for any ipi G £jfe; (i = 1, 2, 3), 

(6.9) < G^i,^),^ >ff=0. 
The direct calculation shows that 

(6.10) G$, V>f Jl011 ) =0 for i> G {^f (23l)00 ,^2 (2Jl)00 >f° 021 }, * = 1,2. 

Then by ^ = a^iOii^f" 011 + 2/7iOiiV , 2 n ° 11 + $(^.71011= 2/ji0ii) and (6.4)-(6.10), we 
derive that 

<G(^),< 31011 > ff 

= < G(V>? ,l0U ,*),*?' 10 " >h ^011+ < G(^ l011 ,*),< Jl011 >h tfjxon +o(3), 
=- < G(vf ,l0U J ^ lG "),* >« ^1011- < G(v£ ,l01 \*?* ou ),* >« feon + o(3), 

2^iGi^ 6 a a A 2 C 2 7r 2 , 2^2, , ,„x 

(^iion - 2/^011)^1011 - „ ~ (gjion + l/jiOiiJ^iiQii + °( 3 )> 



<raia 2 of !tfl)0 v " JlUil ^ 8a lC * 2 

2AiGi7r6 

i ^%'iOii2/iiOiij2/jiOii 

craia 2 a (2ii)0 

/ g^G^ A 2 G 2 7r 2 w 2^2, , ,,x 

= - ( aa 1 a 2 4 Ji)0 + ta^"^*™ + »*0" )**«»• + ° (3) ' 

Similarly, we obtain 

, n( , ,s .T.ftlOll . / 2AlGl7T 6 ^42G 2 7T 2 2 ,2s , ro\ 

< Gfcfc tf), * 2 > H = -( ^^ + -^7)^1011 + 26iOii)2/;iOii + 0(3). 



Hence, the reduction equations are given by 
(6.11) 



dt 

/3jiOii(-R)2/3iOii + S(x 2 ji011 + 2/Ii0ii)%i0U + o(3), 



dt 



^2/jlOll _ , , f>! x , ;i 

where 

(6.12) S = + + AlCl ) + l + A2G2) < 0. 

CTa ( 2jl )o 8 Ji a i 

A standard energy estimate on (|6.1ip together with the center manifold theory show 
that, for R < R Cl , (U, T) — is locally asymptotically stable for the problem (|2.1[) - 
3J. Hence by Theorem fCTl the solutions to bifurcate from (U, T, R) = 
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(0, R Cl ) to an attractor Sr. Moreover, by (|6.11|) - (|6.12|) together with Theorem 5.10 
in [18j . we conclude that Er is homeomorphic to S 1 in H. 

6.2. Completion of the proof of Theorem 3.4. In this subsection, we prove 
that Er consists of steady state solutions. It is clear that the first eigenvalue of 
Lp\jj^ is simple for R rj R Cl . By the Kransnoselski bifurcation theorem (see among 
others Chow and Hale [4] and Nirenberg [21]), when R crosses R C1 , the equations 
bifurcate from the basic solution to a steady state solution in H . Therefore the 
attractor Er contains at least one steady state solution. Secondly, it's easy to 
check that the equations (12.ip ~ (|2.5p defined in H are translation invariant in the 
cc-direction. Hence if ipo(x, y, z) = (U(x, y, z), T(x, y, z)) is a steady state solution, 
then ipo( x + Pt y, z ) are steady state solutions as well. By the periodic condition in 
the x-direction, the set 

SVo = {^o(x + p, y, z)\p s K} 

is a cycle homeomorphic to S 1 in H. Therefore the steady state of (I2.ip - (|2.5j) 
generates a cycle of steady state solutions. Hence the bifurcated attractor 
consists of steady state solutions. The proof of Theorem 13.41 is complete. 

6.3. Proof of Theorem 3.5. The proof follows directly from the classical Hopf 
bifurcation theorem and Lemma 15.61 
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